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Abstract 

Recent experiments on the laser-induced resonant annihilation 
provide a precise measurement of transition energies of antiprotonic 
helium atoms. Relativistic corrections of an order of a 2 to the pure 
Coulomb interaction will be taken into account in the theoretical 
description of energy spectra of antiprotonic helium atoms. The 
splitting of energy levels arising due to the spin-dependent part of 
the relativistic interaction is considered for a number of states of 
3,4 Hepe systems. 
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1 Introduction 



Metastable antiprotonic helium atoms 3,A Hepe have been discovered in ex- 
periments on the delayed annihilation of antiprotons in helium media 
||. Analogous long-lived systems were observed in experiments with nega- 
tive kaons and pions ||. The discussion of the theoretical calculations 
on antiprotonic helium atoms and related topics can be found in ||. 

The precise measurement of transition energies of antiprotonic helium 
atoms in recent experiments on the laser-induced resonant annihilation 0, 
fTj, || invokes the theoretical description of energy spectra with comparable 
accuracy. Such a description of energy spectra requires that minor effect 
of relativistic and QED interactions and the coupling with the continuous 
spectrum should be taken into account. The relativistic corrections of an 
order of a 2 (a = e 2 /Tic - fine structure constant) to the pure Coulomb 
interaction are the most important and should be firstly considered. Next 
in line are QED corrections to energies of higher orders in a. 

Since the contribution to energies from the relativistic interactions de- 
pends on the antiproton mass, charge and magnetic moment, the compari- 
son of precise calculations and measurements of the energy spectra can be 
used for determining the antiproton properties. This knowledge is signifi- 
cant in testing the fundamental symmetry principles. The detailed discus- 
sion of this problem can be found in ||. 

The spin-dependent part of the relativistic interactions gives rise to split- 
ting of energy levels, and each single transition turns into a multiplet. Suf- 
ficiently large distances between lines in the multiplet can be measured 
experimentally. It is worthwhile to mention that the resolution in current 
experiments is about 10GHz and without much difficulty can be improved 
to 1GHz [[!(]]. As it will be discussed below, due to the interaction with 
electron spin, antiprotonic helium energy levels split into two multiplets 
and the interaction with nuclei spins provides a minor splitting into each 
multiplet. Values of the former large splitting are presented in this report. 

The relativistic interaction in antiprotonic helium atoms is described in 
the next section, and in section 3 the method of calculation of the level 
splitting is discussed. Numerical results are presented in section 4; outlook 
and discussion, in the last section. 



2 Relativistic interaction 



For each pair of particles i, j in the three-body system the relativistic cor- 
rection of an order of a 2 to the pure Coulomb two-body potential can be 
described by the Breit interaction of the form 



Uij = -aPZiZjfymi 2 + m i 2 )5(r^) + 2m ^ (PiPj + /-/i^ir/j), )i> ; ) + 
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where = r, — r 3 -, r i; p i; s i; mj, Zj, /Xj are the radius-vector, momentum, 
spin, mass, charge and magnetic moment (in units eh/2mic) of particle z. 
Here and below atomic units (h = e = m e = 1) are used. The correction to 
the kinetic energy of an order of a 2 for each particle i is 
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Full relativistic correction H r of an order of a 2 to the three-body nonrela- 
tivistic Hamiltonian is a sum of Uy for all pairs of particles and ATj for all 
particles 

//,. y^AT, ■ Y^r,,. (3) 

t i>j 

Using expressions U^, ATj in eq. (f|), particles momenta p« will be taken in 



the center of mass frame of the three-body system [11]. Below, for definite- 
ness helium nucleus, antiproton and electron will be enumerated as particles 
1,2,3, respectively. Relative coordinates r = r 2 — ri, p — r 3 — ri, corre- 
sponding momenta p = — zV r , q = — iVp and angular momenta 1 = [rp], 
A = [pq] will be used to simplify notation. 

The interaction H r , given in @, conserves the sum J = L + J2i s i °f the 
total angular momentum L = 1 + A and particle spins Sj. Each level of the 
nonrelativistic Hamiltonian splits into four and eight sublevels for 4 Hepe and 
3 Hepe systems, respectively. Due to very small mass ratios ms/rrii, ms/m2, 
the largest contribution to the energy splitting comes from the interaction 
with the electron spin s 3 . Taking into consideration only terms responsible 



for the splitting in (0), (|3|), this part of relativistic interaction can be written 
as follows: 



H s = a 2 (^As3 + - ([r - p, q]s 3 ) - (4) 
P P\ 
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While the last two terms in (g) are inversely proportional to the heavy 
particle masses m^2, their contribution to the energy splitting is nevertheless 
comparable to the contribution from the first two terms due to the following 
reasons. The small mass factor is compensated in part due to the large 
angular momentum I ~ L of heavy particles. At the same time, only small 
components of the wave function corresponding to the nonzero electron 
angular momenta A ^ lead to the nonzero splitting value from the first 
two terms in (|5f). 



3 Level splitting 

The interaction H s , given in (J|), conserves the sum j = L + S3 of the total 
angular momentum L = 1 + A and electron spin S3 and splits each level 
into two sublevels, corresponding to the eigenvalues j — L± 1/2. The part 
of the interaction depending on heavy particle spins removes the remaining 
degeneracy and splits each j = L ±1/2 sub level further into two or four 
levels for the A Hepe and 3 Hepe systems, respectively. Values of this secondary 
splitting are much smaller in comparison with the splitting, arisen due to 
the interaction with the electron spin By this reason only calculation 
of major splitting will be presented in this report. 

The nonrelativistic Hamiltonian of the antiprotonic helium atom is 
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H = -— A r - — A p V r ■ Vp + (5) 
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where l//z=l/m 1 + l/m 2 , l//ii = l/m 1 + l/m 3 . The nonrelativistic wave 
function ip LN and energy E LN is the iV-th solution and eigen-energy of the 
Schrodinger equation 

{H-E ln )tP ln = (6) 

for the total angular momentum L. Since the splitting is small in comparison 
with energy differences between states of different L values, the energy shift 



Aln can be found in the first order of perturbation theory over H s 



A jLN = (^ jLN \H s \^ lLN ), (7) 

where ^ jLN is the vector production of tpLN and spin function describing 
the dependence of the electron spin. 

Since the interaction H s (|5]) is of the form H s = As 3 , the energy shift 
Aj LN can be expressed 

j(j + l)-L(L + l)-3/4 ,, 

A jLN = ; — Wln\\A-\\v>ln)i (8) 

2JL(L + 1)(2L + 1) 



where the notation (•[[ • [[•) stands for the reduced matrix element. Level 
splitting AEln = Al+i/2LN — Al~i/2LN is a difference of shifts for the 
j=L± 1/2. 

Due to smallness of the relativistic interaction, radiative transitions pro- 
ceed only between states of the same j. By this reason, in experiment each 
spectral line of the transition from the state LjiVj to state LfNf is to be 
split into a doublet with the interline distance Av = AEl i n 1 — AEi f ^ r 



4 Numerical results 

The variational method, described in [12], was applied to determine eigen- 
functions and eigenenergies of the Schrodinger equation The set of 
simple variational trial functions of the form 

Xnmii*, P) = ya M (i, p)r w p x exp(-a n r - b k p), (9) 

where 3^ A/ (f,p) are bispherical harmonics of angular variables, was used 
in the calculations. 

Splitting values AE LN have been calculated as described above (0), (|^) 
by using variational nonrelativistic wave functions. Up to 600 trial func- 
tions (reftrial) containing up to 15 bispherical harmonics were used in these 
calculations. Nonlinear parameters a n , bk were taken the same as in the pre- 



vious variational calculation of energies and radiative transition rates ||12|| . 

Splitting values for a number of states of the 3 ' 4 Hepe systems in the 
range of experimentally observed values of the total angular momentum L 
are presented in Table 1. 



Table 1. Splitting values AEln (10 6 au) of the lowest levels in the 
3 ' 4 Hepe systems. 



4 Hepe 


N 


L=32 


L=33 


L=34 


L=35 


L=36 


L=37 


1 


-1.10 


-1.15 


-1.15 


-1.14 


-1.12 


-1.09 


2 


-1.12 


-1.09 


-1.08 


-1.07 


-1.04 


-1.00 


3 


-1.01 


-1.02 


-1.00 


-0.98 


-0.94 


-0.90 


4 




-0.94 


-0.94 


-0.90 


-0.86 


-0.82 


5 






-0.93 


-0.90 


-0.84 


-0.81 


3 Hepe 


N 


L=31 


L=32 


L=33 


L=34 


L=35 


L=36 


1 


-1.20 


-1.16 


-1.19 


-1.19 


-1.18 


-1.14 


2 


-1.14 


-1.19 


-1.15 


-1.12 


-1.08 


-1.04 


3 


-1.08 


-1.06 


-1.05 


-1.04 


-1.00 


-0.98 


4 




-0.97 


-0.92 


-0.86 


-0.85 


-0.81 



As it follows from expression (^j) , the form of the wave function at small 
interparticle distances is the most important in evaluating the integral (|7|). 
Convergence of the calculated splitting provides a few per cent relative 
accuracy. It is worthwile to mention that due to the variational method 
of calculation the accuracy is better for the large L and small N states. 
It is impossible to trace the convergence in the case of short lived states 
due to a small multipolarity Al < 3 of the Auger decay. This problem is 
closely connected with a large natural width of these states which exceeds 
significantly a splitting value. Also, the variational procedure meets some 
difficulties in describing the short range behavior of the wave function for 
the large enough N states, especially in the 3 Hepe system. These are reasons 
to omit the above-mentioned cases in Table 1. 

The last two terms in eq. (|5|) describe the interaction of the electron 
magnetic moment with the magnetic field of heavy particles. These terms 
give rise to the largest contribution to the energy-level splitting. For bet- 
ter understanding the splitting dependence on L, N this contribution is 
presented in Table 2 for the 4 Hepe system. The contribution to the energy- 
level splitting from the first two terms in (^) are of the opposite sign and 
much smaller in magnitude. Nevertheless, decreasing in this contribution 
with increasing L compensates the L dependence of the last two terms in 
eq. (|5[) and provides a very slow dependence of the total splitting AEln on 
L. 



Table 2. Contribution of the last two terms in eq. (|5|) to the energy-level 
splitting AEln (10 _6 au) in the A Hepe system. 



N 


L=32 


L=33 


L=34 


L=35 


L=36 


L=37 


1 


-1.41 


-1.43 


-1.40 


-1.37 


-1.34 


-1.28 


2 


-1.39 


-1.34 


-1.30 


-1.27 


-1.22 


-1.16 


3 


-1.26 


-1.24 


-1.20 


-1.15 


-1.10 


-1.04 


4 




-1.14 


-1.11 


-1.06 


-1.00 


-0.94 


5 






-1.10 


-1.06 


-0.98 


-0.95 



5 Discussion 

Due to almost exact conservation of the j value in the radiative transition 
the spectral line splitting will be found as a difference of AEln presented 
in Table 1. Most appropriate for the experimental measurement are the 
favoured transitions between states of the same N, which have the largest 
radiative rates 12]|-|l4"|. However, the calculated splitting values are almost 



independent of L for a given N, and it is not plausible to resolve such a 
small difference in splitting for the favoured transitions. For this reason, 
the experimental proposal for the near future |IJ is aimed at searching for 
the splitting in unfavoured transitions (L, N) — > (L — 1, N + 2). 

In order to measure splitting in experiments on the laser-induced reso- 
nant annihilation the initial state will be long-lived. This is provided by the 
condition that the multipolarity of the Auger decay for this state is Al = 4. 
The next condition is that the natural width of the short-lived final state 
will be smaller than the splitting value, and the multipolarity of the Auger 
decay for this state will be Al = 3. The spectral line splitting for a number 
of suitable transitions is presented in Table 3. These values are of an order 
of the experimentally measurable value ~lGHz. 

Comparing the splitting values for the A Hepe and 3 Hepe systems one can 
mention in the A Hepe case a slower decreasing in AE^ with increasing N. 
As it is clear from Table 3, this isotopic effect is also conserved for the 
spectral line splitting Au. 

Table 3. Spectral line splitting Az/ = AE^n* — AEl } n } (GHz) for the 
transitions E L . N . — > E LfNf in the 3 ' 4 Hepe systems. 



4 Hepe 


3 Hepe 


LiNi->L f N f Av 


LiN^LfNf Av 


33,1^32,3 -0.92 


32,1->31,3 -0.53 


34,1-^33,3 -0.86 


33,1^32,3 -0.86 


34,2^33,4 -0.91 


33,2^32,4 -1.22 


35,2^34,4 -0.88 


34,2^33,4 -1.35 


35,3-^34,5 -0.34 





The following considerations can be used to understand qualitatively the 
L, iV-dependence of the energy-level splitting. Contribution to splitting 
from the interaction of the electron magnetic moment with the magnetic 
field of heavy particles is described by the last two terms in the splitting 
interaction H s (|5]). This contribution is proportional to the relative momen- 
tum of heavy particles p. One can consider that the motion of heavy particles 
is approximately the same as in a hydrogen-like atom and momentum p is 
inversely proportional to the angular momentum L. This is the reason for 
increasing this contribution with decreasing L, as presented in Table 2. The 
contribution from the first two terms in the splitting interaction H s is con- 
nected with the electron rotation and proportional to the small component 
of the wave function arising due to polarization of an electron by p. With 
decreasing L the antiproton moves to a region of increasing electron density 
and the polarization increases. In such a way contributions to the energy- 
level splitting from the last two terms in H s and remaining part of splitting 
interaction are of opposite signs and level off the dependence of the total 
splitting AE LN on L. 

One can consider quasiclassically that the antiproton orbit became more 
stretched with increasing N at fixed total angular momentum. By this 
reason all the terms of the splitting interaction H s decrease with encreasing 
L and provide the N dependence presented in Tables 1, 3. 
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